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Abstract 



An overview of recent studies of nonequilibrium bound interfaces is 
given. Attention is focused on Kardar-Parisi-Zhang interfaces in the 
presence of upper and lower walls, interacting via short- and long- 
ranged potentials. A comparison with equilibrium interfaces is carried 
out, and connections with other nonequilbrium systems are illustrated. 
Experimental realizations of the phenomenology described in this ar- 
ticle are briefly discussed. 

1 INTRODUCTION 

An interface is the moving or static boundary between two distinct bulk 
phases. Effective descriptions that focus on the statistical properties of the 
interfacial degrees of freedom have proved particularly useful, and over the 
past decades received a great deal of attention [H 12 El El • Within this ap- 
proach, the essence of the bulk contributions is incorporated into the (con- 
tinuous or discrete) model parameters; other effects, such as geometric con- 
straints, are modeled as nontrivial surface terms. Examples of the latter 
abound, and include the morphology and motion of an interface constrained 
by the presence of a third phase, by an external field, by inhomogeneities of 
the medium, etc. In thermal equilibrium, or near equilibrium, enormous the- 
oretical, as well as experimental, progress was made. For example, wetting 
in the two-dimensional Ising model was shown to be uniquely described by 
an interfacial model on sufficiently large length scales jS]. Of more practical 
interest are three-dimensional systems with realistic interactions, for which 
a number of theoretical predictions have been confirmed experimentally [H]. 
However, for nonequilibrium systems this field of research is only just start- 
ing. 

In this review, we give an overview of recent theoretical results on nonequi- 
librium bound interfaces. The focus will be on a few models, mostly in the 
form of continuum stochastic growth equations, that are expected to be rep- 
resentative of a broad class of physical phenomena. The primary reason for 
this choice is that, the statistical properties being identical to those of the 
discrete models, continuum equations are more adequate to describe the col- 
lective behavior at the macroscopic level, and thus give a unified picture of 
the field. Although the discussion is intended to be general, unless otherwise 
stated, results pertaining to one-dimensional systems only are given. 
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This paper is organized as follows. Firstly, in Section 2 a brief discus- 
sion of the kinetics of two representative, free-interface models is provided. 
In Section 3, we review the implications brought about by the inclusion of 
interacting, binding walls. Next, we explore analogous phenomena in other 
systems and in the last section describe experimental systems where bound- 
KPZ behavior could be observed. Much of this material has been covered 
from a different perspective in a recent review [7] . 

2 TWO MODELS OF FREE INTERFACES 

2.1 The Edwards- Wilkinson equation 

The following Langevin equation, referred to as the Edwards- Wilkinson equa- 
tion, is commonly employed to describe on a coarse-grained level the dynam- 
ics of an interface separating two bulk fluid phases 

^^ = z.V^/^ + r^(x,t), (1) 

where /i(x, t) is the local height of the interface above some reference plane. 
u is the interfacial tension or the interfacial stiffness if any of the coexisting 
phases is anisotropic. It measures the strength of the Laplacian, a relaxing 
term that arises from bulk evaporation-condensation processes. Microscopic 
fluctuations, which in the present case are assumed to have a thermal origin, 
are modeled by the stochastic noise term (77) , with mean (77) = and variance 
{ri{-x,t)ri{-x',t')) = 2D6(t — t')6{'x — x'), where (■) denotes an average over the 
noise distribution. Since the system is in thermal equilibrium, the noise 
correlation function is related to the temperature through the fluctuation- 
dissipation theorem, i.e. D = ksT. Finally, note that this model ignores 
configurations containing overhangs and bubbles, as /;.(x, t) is a single- valued 
function. 

An interface governed by the EW equation does not move on average 
irrespective of its initial position. This agrees with the thermodynamic pic- 
ture that at bulk coexistence any arbitrary fraction of the system may be 
in one phase, with the remainder in the other. Adding a constant /i to (jT]) 
drives the system out of coexistence. For instance, if /i > the steady-state 
interface moves upwards with an average velocity v = ^, thereby favoring 
the y < /i(x) phase over the unstable y > h{-x) one. In this context, fj, may 
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be interpreted as an external force acting on the interface or as the chemical 
potential difference between the two phases. 

The fluctuations of the interface around its equilibrium position are char- 
acterized by the interfacial width, a conventional measure of interfacial rough- 
ness, defined as 

W'{L,t) = {[hit)-h{t)r), (2) 

where h{t) is the spatial average of h{'x, t). An initially fiat surface is essen- 
tially confined between h{t) — W{t) and h{t) + W{t), at later times. It is 
well-known that the behavior of W depends on the system dimensionality: 
for d > 2 the interface is said to be smooth as W converges in the limit 
L,t ^ oo. At d = 2 the interface is logarithmically rough, 

^ lln(L/a), 

Here a is a microscopic scale of the order of the bulk correlation length and 
defines a time scale characterized by the so-called dynamic exponent z. 
Finally, in c? < 2 the interfacial behavior is universal and given by. 

The initial increase of the interfacial width as a power of time is governed 
by the growth exponent, f3; this is followed by a saturation regime where 
the interfacial width scales as W{L, oo) ~ L''. The crossover between the 
initial and the saturation regimes occurs when the lateral correlation length 
^{t) ~ t^/^ is of the order of the size of the system. The roughness exponent 
is a measure of the typical transverse fluctuation W of a. given interfacial 
segment L as given by This behavior is summarized by the Family- 
Vicseck scaling relation W{L,t) ~ L''f{t/L^), where the scaling function / 
satisfies f{u) ~ for m ^ 1 and f{u) = const for m ^ 1 Clearly, 
continuity at t = requires ( = (3z. Note that in an infinite system the 
width does not saturate, implying that the interface is rough in the sense that 
it makes arbitrarily large excursions from its average position. Due to the 
linear character of equation (0) the exponents can be calculated in arbitrary 
dimension. In particular, z = 2, [3 = {2 — (i)/4, and C = (2 — d)l2. 

Equation (Q) can be derived by differentiating a simple Hamiltonian 

H{h) = / dx -{Vhf (5) 
Jo i2 J 
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that is recognized as the square-gradient approximation to a surface free 
energy that is proportional to the total interfacial area, 

1/2 

(6) 

The study of systems off-coexistence requires the inclusion of an additional 
term /i / /i(x)(ix in where /i is the chemical potential difference between 
the two phases. Of course, all the stationary properties of the EW interface 
described above can be recovered from (0) in a straightforward manner. 



{Vhf 



2.2 The Kardar-Parisi- Zhang equation 

As referred previously, the velocity of an interface described by the EW equa- 
tion is f = /i, regardless of its initial position. At /i = this is equivalent to 
bulk phase coexistence. A natural generalization of (0) consists in assuming 
that V is no longer constant, but depends on the instantaneous local-slope 
of the interface at x. Thus, substituting /x by viVh) and using again the 
small-slope approximation one finds^ 

v{Vh) ^ t;(0) + Vv{0) ■ Vh + ^[{Vh) ■ V]h{Vh = 0). (7) 

Assuming isotropy in the sense that d'^v / dx-idxj oc and removing the first 
two terms through a Galilean transformation h{'x,t) /i(x + Vv{0)t,t) + 
v{0)t, yields the celebrated Kardar-Parisi-Zhang (KPZ) equation [TTj . 

= uV'h + XiVhY + r/(x, t). (8) 

The only difference between the KPZ and the EW equations is the nonlinear 
term A(V/i)^, that is the most relevant nonequilibrium perturbation to the 
equilibrium EW equation and is expected to describe the scaling properties 
of rough surfaces growing in the absence of conservation laws. Somewhat sur- 
prisingly, experimental realizations of the KPZ behavior are far from evident 
in marked contrast with its presence in a variety of theoretical models. It 
has been argued that this may be due to the occurrence of medium-induced 



^It can be shown that third and higher order terms may be neglected in the large-scale 
limit [ini. 
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non-local interactions ^21, and/or long instability-induced transients that 
mask the asymptotic KPZ scaling [T3]. 

The Family- Vicsek scaling picture (j3)) of the EW equation also holds 
for the KPZ. It turns out that in one dimension the stationary probability 
distribution 

P(h)=exp{-^liVhrdx), (9) 

which is a solution of the Fokker-Planck equation of the linear theory (A = 0) 
in any dimension, is also a solution of the Fokker-Planck equation corre- 
sponding to the KPZ. Hence, the nonlinear ity does not affect the steady 
state solution and consequently the two interfaces have the same roughness 
exponent ( = 1/2. The other two exponents may be obtained from the re- 
lation z = 2 — ( valid in any dimension, and from (3 = a/z, previously 
derived. In dimensions d > 1, the scaling exponents cannot be obtained 
analytically, nor perturbatively. At d > dc = 2, both a weak and a strong- 
coupling regimes occur depending on whether the value X^D/u^ is smaller or 
larger than a critical value. In the weak-coupling regime, which is unstable 
at the critical dimension dc = 2, X vanishes asymptotically leading to the 
EW behavior z = 2,( = 0- By contrast, the strong-coupling regime is con- 
trolled by a fixed point that is not reached by perturbation analysis. Several 
conjectures for the critical exponents have been put forward [T], including a 
recent claim of a mathematical proof j^. Another controversial aspect of 
the KPZ behavior concerns the existence of the upper critical dimension. For 
a more elaborated discussion of these (rather technical) matters, the reader 
is referred to existing reviews El . 

The KPZ nonlinearity is related to the lateral growth that occurs when 
a depositing particle sticks to the first particle it encounters on the surface 
pp. As a consequence, by contrast to EW-like surfaces, interfaces governed 
by the KPZ equation have a nonzero velocity at = 0, 

V = [^{dth)dx = X l^{{yhf)d'^x = Am^, (10) 
Jo Jo 

where m? is the overall quadratic slope of the interface. Since the stationary 
probability distribution is known in d = 1, = —d^lnZ can be computed 
through the simple Gaussian path integral, 

Z = / P/i exp ( - 7 f {Whfdx) . (11) 
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with 7 = v/2D. Thus, v = DX/Aua, which depends on the microscopic 
scale a and is therefore nonuniversal; however, it has a universal finite-size 
correction v{L) = DX/2h'L, a result that holds for t ^ . See [T7] for a 
different, more detailed derivation. 

As the nonlinear term of the KPZ is purely kinetic in origin, the KPZ 
equation cannot be obtained by differentiation of an equilibrium Hamiltonian 
such as (0).^ Nevertheless, it can be mapped onto a linear diffusion-equation 
by a simple change of variables, the so-called Cole-Hopf transform /i(x, t) = 
{y/\) lnn(x, t), yielding 

MZ^ = ^V2n + nr/. (12) 
ot 

This constitutes the directed polymer representation of the KPZ P] . We shall 
take advantage of this representation in a later section but, for now, simply 
point out the presence of the multiplicative noise term. 



3 BOUND INTERFACES 

Let us assume that the interfacial fiuctuations are restricted by the presence 
of another interface, by a physical barrier, or by any other mechanism capable 
of confining it. We will see shortly how to incorporate this effect into a 
growth equation. Consider a rigid, impenetrable wall that prevents large 
interfacial excursions in a given direction. In this case, the interfacial width 
approaches a constant value, W{L,t) as L,t — >• oo even if the free 

interface is rough. Therefore, a bound interface is never rough. Moreover, 
since the fluctuations are cutoff in a given direction, there will be an effective 
fluctuation-induced repulsion between the wall and the interface, leading to 
the possibility of interfacial pinning-depinning (or unbinding) transitions. 

In the following we briefly review the behavior of a modifled EW equa- 
tion describing an interface bound to a wall, for which the aforementioned 
behavior can be worked out explicitly within the mean-fleld approximation 
|19j . The case of the bound KPZ interface will be tackled afterwards. 

^It has been argued, however, that equation © is not correct unless dth and rj are 
projected along the normal direction, in which case a non-bound Hamiltonian exists |18| . 
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3.1 A bound Edwards- Wilkinson equ-ation 

Restricting the height variables in the EW equation to, for instance, positive 
values requires supplementing with the hard wall condition V{h < 0) = 
oo. This is not only inconvenient for analytical studies, but also unnecessarily 
limited in scope. Therefore, more general types of walls will be considered, 
the only requirements being the soft wall condition V{h) —>■ oo as h —>■ h^aii, 
and V{h) — as /i — oo with V{h) differentiable. In this case, the wall 
interacts with the interface, in addition to cutting off some of its fluctuations. 
The large-scale interfacial behavior as well as the character of the unbinding 
transitions are not expected to change as the hard-wall is replaced by a soft- 
wall ^19j. We therefore consider the EW growth equation in the presence of 
an effective interface potential V{h) fIU\ OT] 

— ^ = ^^'h + ^(x, t), (13) 

where rj is Gaussian white noise with the same mean and variance as in (^J. 
The region y > h{x.) corresponds to the stable thermodynamic phase, say 
A, the region y < /;,(x) corresponds to a second phase B that coexists with 
A when /i = and /i(x, t) is the local height of the AB interface measured 
from the wall (see figure (1)). 

ly is the interfacial tension of the AB interface and V{h) accounts for 
the net interaction between the interfaces binding the B layer (the wetting 
film), in this case the substrate and the AB interface. In the simple case 
where AB corresponds to a vapor-liquid interface, a repulsive V{h) simulates 
adsorption of the liquid phase from the vapor. The form of V{h) depends on 
the nature of the forces between the particles in the fluid and with the wall, 
and its derivation from bulk, microscopic Hamiltonians is far from trivial. 
Ideally one constrains the interface, away from its equilibrium fiat position, 
in the configuration h{x.) and using the microscopic Hamiltonian, takes a 
partial trace over the bulk variables. Of course this cannot be done exactly 
for realistic Hamiltonians. In approximate derivations one usually requires 
consistency at the mean- field level. If all the microscopic interactions are 
short-ranged, one may take for sufficiently large h at bulk coexistence, 

V{h) = b{T)e~^' + ce-2\ (14) 

where h is measured in units of the i?-phase bulk correlation length, T is the 
temperature, b{T) vanishes linearly as T — T^;, the wetting temperature (see 
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below) and c > 0. If, instead, one considers long-ranged (van der Waals) 
interactions, the potential has the general form [221 



and in this case b{T) is related to the Hamaker constant. In all cases, if we 
add the contribution from the chemical potential difference /i/i to V{h), the 
behavior of the system will be controlled by the latter term, since this is a 
bulk contribution that dominates over the surface term V{h) except at bulk 
coexistence (figure (2)). Finally, note that the present description assumes 
the existence of an interface and thus it is only valid below the bulk critical 
temperature, where distinct "liquid" and "gas" phases are defined. 

Some remarks concerning the connection of equation (fT^ with wetting 
phenomena follow. Wetting occurs when the wall preferentially adsorbs one 
of the phases, say B, while the bulk may be in a different thermodynamic 
state. At a wetting transition, the thickness (h) of the adsorbed B layer di- 
verges. This occurs at all temperatures above a certain wetting temperature, 
T^, at bulk phase coexistence, i.e. ai fi = fic = 0; by contrast, at /i 7^ /ic 
the thickness of the liquid film may be large, but cannot diverge (bound 
interface). 

Equation (jl3j) is a djTiamic model for the relaxation of the interfacial 
height /;,(x, t) towards the equilibrium configuration that minimizes the Hamil- 
tonian 



in the limit t ^ 00. Within mean- field approximation (h) follows from the 
condition dV{h)/dh = 0, whereby one finds that the equilibrium thickness 
of the wetting layer for an attractive wall {b < 0) and long-ranged forces 
with m = 3,n = 2, (non-retarded van der Waals forces [21]) is given by 
(h) = —2b/3c. A critical wetting transition takes place as 6 ~ T — — > 
bw — 1 i-G- the wetting layer thickness diverges as {h) ~ |6 - b^f- with 
(3c = —1- This mean- field result was observed experimentally |2S] since the 
upper critical dimension of the long-ranged system is less than three [^ . 

Critical wetting driven by short-range forces may also occur since the 
long-ranged interactions can be neglected if the bulk correlation length is 
sufficiently large i.e., when is close to the bulk Tc. This has been con- 
firmed recently by the observation of effective short-range critical wetting at 
the liquid- vapor interface of methanol-alkane mixtures [2Zj. Within mean- 
field theory, {h) = ln(— 2c/6) and, consequently, /3c = 0{log). However, the 



V{h) = b{T)h-"' + ch 



—n 



n > m > 0, 



(15) 




(16) 
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measured exponents, which are also mean-field like, are at odds with the 
theoretical, renormalization-group-based predictions |27j . 

By contrast, the wetting transition may be driven by the chemical poten- 
tial difference between the A and B phases, at any temperature above T^; 
this is always a continuous transition and it is known as complete wetting 
|23j . A study of complete wetting transitions requires adding a linear term 
IJ,h to the Hamiltonian, where n is the chemical potential difference of phases 
A and B. Thus, on approaching coexistence for T > T^^, (6 > at the mean- 
field level), {h) diverges as (h) ~ \fi — /xd"^'', with /ic = 0, and Ph = 0{log) 
and Ph = l/{m + 1) for short- and long-ranged forces, respectively. Exper- 
imental observations of complete wetting transitions are numerous and are 
characterized, in general, by long-range mean-field exponents (at least far 
from the bulk critical temperature). The interested reader is referred to [221 
for reviews, and to and |28j for a detailed account of recent experimental 
and Monte Carlo results, respectively. 

This rich behavior is summarized in the schematic phase diagram of figure 
(3): (along path 1) {h) diverges continuously as coexistence is approached 
from the gas phase (complete wetting); (path 2) as T approaches T^j at 
coexistence (h) may diverge either discontinuously (first-order wetting ) or 
continuously (critical wetting); (path 3) the wall remains non wet when coex- 
istence is reached . Clearly, in the wetting regime at long times, the form of 
V{h), determines the rate of growth that is logarithmic, i.e. (h) ~ t^'^ ~ Int 
OT 9h = for short-ranged interactions and a power law, (h) ~ ti/(2+»")^ for 
long-ranged repulsive interactions decaying as /i"™ Off coexistence, for 
a small and negative 6fi = fi — fic, the wetting layer thickness grows with 
time as predicted by the mean-field theory, provided that t is much smaller 
than the correlation time 6fi'^^. Then, it relaxes exponentially to its equilib- 
rium value {h{t = oo)). These exponents are related by 9h = Ph/zi'x, and 
Ut = 120] • On the other hand, if 5/i is positive, {h{t)) grows linearly. 

The results for the exponents are summarized in table (0) and Monte Carlo 
simulations of the growth of wetting layers can be found in [211 ISOl • 

The behavior below deserves some comments. Wetting only occurs at 
coexistence. The first-order transition that takes place at /i = has to be 
distinguished from wetting because it is not driven by the substrate. In fact, 
the wall is not wet at the transition point and the depinning transition that 
occurs upon crossing the boundary line /i = is, as such, trivial since to the 
right of the fi = line liquid is the only stable thermodynamic phase, with 
or without a wall. 
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Lastly, we comment briefly on how the effective wall interaction due to 
fluctuations comes about. As a result of the conflnement of the interface, 
there is an increase in the elastic bending energy and an entropy loss. For 
d = 1 both contributions are of the same order and can be estimated as 
Ae ~ TAs = AV ~ while for d = 2 one has exp(— 2/i^) ^H]- By simply 
adding AV to V{h) one gets a description of the effects of fluctuations within 
mean-fleld theory. 



3.2 A bound Kardar-Parisi- Zhang e-quation 

Consider the following dynamic interfacial model 

The fact that the KPZ equation is not invariant under the transformation 
h — >■ —h implies that the sign of A assigns an orientation to the wall, distin- 
guishing two types: upper and lower walls. An upper (lower) wall restricts 
(or even prevents) large interfacial excursions into the region h > {h < 0). 
Thus, for negative values of A the interface is pushed on average against a 
lower wall, while for positive A it is pulled away from it, and exactly the op- 
posite occurs at an upper wall. As will be shown below, for a given sign of the 
non-linearity upper and lower walls lead to quite different phenomenologies, 
but the case A > and a lower wall is completely equivalent to A < and an 
upper wall [SI]. We stress that, were it not for the presence of the nonlin- 
ear term, such distinctions would not have been necessary (both equilibrium 
wetting and dewetting are symmetric phenomena). 

As noted earlier, in the absence of the limiting wall, bulk coexistence no 
longer obtains at /i = 0. Rather, a nonzero chemical potential yUc ~ A is 
required to balance the force exerted by the nonlinear term on the tilted 
regions of the interface. For A = the model reduces to the equilibrium one 
and /ic = as usual. As in the previous section, we distinguish the cases of 
short- and long-ranged interactions. Also, for deflniteness, in what follows 
the wall will always be assumed to be a lower one, while A may take either 
sign. 



3.2.1 Short-ranged forces 

Consider the Langevin equation ()17p where V{h) = he~^^ + ce'"^^^. The 
parameter p controls the nature of the wall: the limit p —>■ oo corresponds 
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to an impenetrable wall and the sign of p determines whether it is an upper 
(-) or a lower one (+). The same equation with c = 0, 6 > (repulsive 
wall), and A < was first studied by Tu et al. [32^, and by Muiioz and Hwa 
[31j for arbitrary A in the context of nonlinear diffusion with multiplicative 
noise (see later in this article). The effect of inert and attractive walls was 
then investigated by Hinrichsen et al. in 33, 34j using a discrete solid-on- 
solid model with dynamics that violate detailed balance. The model includes 
particle adsorption and desorption rates, plus an additional growth rate at the 
wall that simulates a short-ranged attractive or repulsive interaction between 
the wall and the interface. Their results were subsequently confirmed and 
expanded by Giada and Marsih who mapped (fT7|) to a Langevin equation 
with multiplicative noise and analyzed it within mean-field, and by ourselves 
jSni EZI, by means of direct integration of (fTTjl with A < 0. The work of [33] 
privileges the role of the noise as the driving force, while in |36| the noise 
strength is fixed and b and are taken as control parameters. The case A > 
was recently investigated in pHj . 

Some of the exponents defined in table (0) niay be related to KPZ ex- 
ponents. In particular, the dynamic exponent z is unchanged from its KPZ 
value, and Vx = 1/(22 — 2) for any value of A [23E2I- Thus, in one dimension 

= 1 since z = 3/2. The two remaining exponents, j3h and 9h are related 
by the scaling form 9h = Ph/^xZ- 

The long-ranged repulsion exerted by the wall on the interface, mentioned 
previously, may be estimated using a simple scaling argument: the wall makes 
itself felt when the characteristic distance of the interfacial excursions is of 
order h, C,± ~ h. From Q and the definition of u^, we find ~ ~ 5^'^"^^, 
leading to an effective repulsive force /i"^/*^^^ ~ /i"^/'^"', since C, = l/2m.d = 1 
for both KPZ and EW. By substituting = 2/3 (see table (0)) we recover 
the EW resuh, while z/^ = 1 yields V{h) ~ h'^ for the KPZ IMj- 

A < 0. We first consider a negative KPZ non-linearity. Before discussing 
the results, we note that upon inverting the sign of h, A changes to —A and 
p to —p, so that the exponential now acts as an upper wall and A > pushes 
the interface against it. This shows that the cases A > — upper- wall and 
A < — lower- wall are equivalent. 

The associated phase diagram in the b — fi plane is depicted in figure (5). 
The solid line is the continuous phase boundary between depinned and pinned 
phases. Between the dashed lines, both of which correspond to first-order 
boundaries, the pinned and depinned phases coexist as stationary solutions of 
the dynamical equation. The three lines meet at the tricritical point {fic, b^)- 
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The unbinding transition at /i = /ic for any 6 > 6^ is the analogue of 
the equihbrium complete wetting phase transition (path 1 of figure (4)). 
The exponent governing the divergence of {h) is (3h ~ 0.41, with error bars 
that exclude the equilibrium value (3h = 1/3 |H3lini- Along path 2, which 
corresponds to a nonequilibrium critical wetting transition, h is progressively 
increased until the nonwet phase becomes unstable at 6^. It is found that 
the critical temperature is depressed from its mean-field value = 0, with 
an associated critical exponent (3c ~ —1.2 57. 

As for the behavior below (attractive wall) there are two first-order 
dep inning transitions by varying fi (paths 3 and 4 of figure (4)). Recall that, 
if the wetting temperature is to be a meaningful concept, the depinning 
transition for b < cannot be viewed as a wetting transition or, in other 
words, these transitions occur not on approaching coexistence, but when 
crossing the boundary line. For sufficiently negative values of b the pinned 
phase becomes unstable at /i*(&) (path 4), and the depinned phase at /ic < 
/i*(6) (path 3). Consequently, in the range /ic < < /i*(&) both phases 
coexist in the sense that if the interface is initially close to the wall it remains 
pinned, while if it is initially far from the wall it detaches and moves away at 
a constant velocity. As a consequence of the "broad" phase boundary there 
is the possibility of defining critical wetting along a range of different paths, 
delimited by the dashed lines of figure (4). It has been checked that the 
exponent (3c does not change when the tricritical point is approached along 
different paths within this region |37j . 

The fact that the coexistence region is finite rather than a line is a 
nonequilibrium effect since at equilibrium two-phase coexistence only oc- 
curs on a surface of dimension one less than the dimension of the space of 
thermodynamic parameters. An equilibrium system initially in a state other 
than the equilibrium one will be destabilized by local fluctuations in the form 
of droplets of the equilibrium phase. To ensure generic multistability, a ro- 
bust mechanism for eliminating these droplets of the minority phase must 
exist. In the present case, such a mechanism was found by Hinrichsen et 
al. [33]: when, due to thermal fluctuations, a segment of the interface over- 
comes the potential barrier and gets out of the potential well (see figure (2)), 
it rapidly acquires a triangular shape, after which it shrinks at a constant 
velocity driven by the negative nonlinear term, in a time proportional to its 
size. Typical interface profiles resulting from numerical solutions of (fT7j) in 
d = 1 and 2 are shown in figure (5) (similar configurations are described in 
the discrete model of reference [33 )• The largest size of the depinned regions. 
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i.e. the size of the triangular bases, increases as the instabihty threshold /i* 
is approached. Once the last site has detached, the interface "takes off" and 
a transition from a pinned to a moving phase takes place. In fact, owing to 
finite-size effects the only stable phase within the coexistence region is the 
moving one, making the analysis of the pinned regime rather hard. This diffi- 
culty is circumvented by studying r, the time taken by the interface to depin 
in the limit L oo. For fi > fi*, t saturates with increasing system-size 
and thus the interface detaches in a finite time. For /ic < < /i*, however, r 
grows exponentially with the system-size so that the pinned phase becomes 
stable in the thermodynamic limit. The stationary distribution of triangles as 
a function of their base-size turns out to be exponential. This indicates that 
there is a typical size for the depinned regions and rules out the possibility 
that growth is driven by a coarsening mechanism ^7\. Another interesting 
feature is that, for a given pair of parameters /i and b, all the triangular 
facets have the same slope, s, which can be determined through the relation 
|Aij|s^ = yU, where is the renormalized nonlinear coefficient of the KPZ. 
Historically, the possibility of a finite phase boundary was first discussed in 
Toom's north-east-center voter model Other examples include systems 
of harmonically-coupled, identical nonlinear constituents under the simul- 
taneous influence of additive and multiplicative noise |12], and a Leshhorn 
automaton subject to transient velocity-dependent forces jlSl. The relation 
between these systems and condiditons for generic phase coexistence are dis- 
cussed in jH]. 

It has been pointed out that upon diminishing both the depth of the 
potential well and the surface tension, the first-order phase transition at n* 
becomes continuous with directed percolation (DP) critical exponents |45j.^ 
This is in line with a recent claim that first-order phase transitions in one- 
dimensional nonequilibrium systems with fluctuating ordered phases are im- 
possible, provided there are no conservation laws, long-ranged interactions, 
macroscopic currents, or special boundary conditions jH^. If this is the case, 
the reported first-order behavior will be a very long transient effect. The 
presence of DP behavior is now a well-established result that was found in a 
variety of discrete models IT7| EH] • Nevertheless, we cannot discard the 
presence of a tricritical point for sufficiently negative values of b, separating 
DP from first-order transitions. At present, this question remains open. 

■'To observe DP behavior a different order-parameter, namely the density of pinned 
sites, has to be used. More on this later 
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A > 0. Consider A > 0. Early results were reported in [SI], but the 
observed critical exponents are far from their true asymptotic values. Later 
studies revealed marked differences when compared to the system with a 
negative nonlinearity and attractive walls IHH] . These are apparent in the 
phase diagram of figure (4), similar to that of A = (figure (3)). Before we 
proceed, we note that the results of this section were obtained using Monte 
Carlo simulations of discrete models. Attempts at numerical integration of 
a bound KPZ equation failed to reach a stationary state. It is well known 
that the results from numerical integration may disagree with the predictions 
from the continuum KPZ ^^1- The usual way to handle the integration of 
the KPZ equation with short-ranged forces and negative nonlinearity, the 
Cole-Hopf transformation, is also plagued with instabilities. A mean-field 
analysis, however, exits for the oder-parameter n = exp{—h) (see section 4). 

We start by describing the results for a repulsive wall. A continuous 
wetting transition is found as yU — > /Xc when h > (6^ = at the mean- 
field level; path 1 of figure (4)). The behavior of the order parameter in the 
vicinity of the critical point yields [3h = 0.52(2), and the time evolution of 
{h) at behaves as {h{t)) ~ t^'^ , with 6h = 0.355(15) j2H]- From the scaling 
relation z = Ph/j^xOh, z = 0.52(2)/0. 355(15) = 1.5(1) in good agreement 
with the theoretical prediction z = 3/2. These results, however, disagree 
with those of reference ^3] probably due to the presence of extremely long 
transients in the latter. The results are collected in table (2). 

For attractive walls, we note that there is a line of first-order depinning 
transitions ending at a tricritical point at b^] (path 3 of figure (4)). By 
contrast to what happens for A < 0, phase-coexistence is restricted to a 
line since the mechanism responsible for the elimination of droplets of the 
depinned phase works only for negative A. Finally, the tricritical behavior 
associated with the transition along path 2 has not been investigated yet. 

3.2.2 Long-ranged forces 

Next, consider the KPZ equation in the presence of a long-ranged potential 
V{h) = b{T)h~"^ + ch~"'. The resulting equation is not amenable to analyti- 
cal treatment and therefore we are limited to numerical results from Monte 
Carlo simulations of one- dimensional, discrete models. Recently, Lipowski 
and Droz |17j studied an interfacial growth model with dynamics inspired by 
the synchronization transition in coupled map lattices and included a power- 
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law interaction between the interface and the wall.'^ The model, however, 
lacks a surface tension and its microscopic rules are such that it is difficult 
to ascertain whether it may be described by a bound KPZ equation. 

Interestingly enough, the KPZ nonlinearity appears to be irrelevant above 
the wetting temperature bui for any sign of A, and an equihbrium complete- 
wetting transition is found along paths 1 of figure (4) j^. The associated 
exponents are collected in table (0) and have been verified numerically in 
a simulation of a variant of a model introduced in |Tn|, whose continuum 
counterpart is known to be the KPZ with A < and a lower wall [5T]. 

Below by^, a phase diagram similar to that for short-ranged interactions 
is found (figure (4)) for either sign of A. In particular, for A < and for the 
set of parameters explored, the transition at (path 4) is characterized 

by DP exponents when the order parameter n = exp{—h) is used For 
A > 0, to our knowledge, there are no pubhshed results. Results of ours 
indicate a phase diagram similar to the one shown on the right of figure (4): 
along path 3 the system undergoes a discontinuous transition. The analogue 
to a nonequilibrium critical/first-oder wetting transition along path 2 has 
not been investigated. 



4 RELATED PHENOMENA 

In this section, the directed-polymer representation of the KPZ is exploited 
to analyze bound KPZ interfaces. As shown in section 2.2, the KPZ can be 
linearized by means of the Cole-Hopf transformation, n(x, t) = exp[a/i(x, t)], 
at the cost of introducing a multiplicative-noise (MN) term. This transforms 
the KPZ equation (jH)) into 

= uV^n + nri. (12) 

at 

Note that the noise amplitude vanishes for n = 0, and neither the determin- 
istic dynamics, nor fiuctuations can take the system out of this state. For 
this reason the configuration n = is known in the MN literature as an ab- 
sorbing state, while n 7^ is commonly referred to as the active phase. This 
representation has been successfully applied to the case of short-range inter- 
actions with a negative non-linearity, enabling analytic calculations both at 

^The relationship between synchronization transitions and a bound KPZ wiU be seen 
in the next section. 
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the mean- field level \6b\ \'6b\ n7\ I42j. and beyond via a renormalization group 
approach [321 EH] • Let us consider again (fTTj) in the presence of a short-ranged 
potential with, for simplicity, —A = i/ = D = 1 (different coefficients for the 
Laplacian and the KPZ nonlinearity may be accounted for by a proper choice 
of a). The change of variables h = — Inn leads to 

9ra(x, t) dV(n) 

with V{n) = fin'^/2 + + and where we have made use of the 

Stratonovich calculus 1^.^ In general, potentials of the form V{n) = 
jjj^p+'i _j_ cfi'^p+'^ with p > result in equivalent effective Hamiltonians since, 
when expressed in terms of h, p can be eliminated by defining the length 
scale. In the n language, unbinding from the wall, (h) oo, corresponds to 
a transition to an absorbing state, (n) 0, recovering the usual behavior of 
Landau theory where the order parameter becomes small rather than large 
near the transition.^ To complete the analogy between the MN and the 
interfacial descriptions, one can give a physical interpretation to n by noting 
that it is the density of sites at zero height, n{'x,t) = Sh(:x.,t),o The 
behavior of (n) is characterized by a set of exponents similar to that of (h): 
(n) ~ t~^", (n) ~ ^ ~ Sfx~'^'^, and ^ ~ t^^^. Some of these exponents are 
obtained from the KPZ exponents using scaling arguments. For instance, it 
was shown in [321 EH] that the dynamic exponent for the MN case is identical 
to the value of the KPZ z, that is related to it through = l/{2z — 2), and 
that > 1. Also, the scaling relation 6n = (3n/i'xZ is satisfied. It is therefore 
not surprising that the phase diagram of the MN equation (jl8p is very similar 
to that of the KPZ. There is a strong noise phase for d < dc = 2, and both 
a weak and a strong noise phases for o? > 2 depending on the noise intensity 
D. Weak and strong phases obtain for values of D that are, respectively, 
smaller and larger than a critical value |22]- See j^S] for a comparative list 
of critical indices below and above dc- 

In this context, (jl8|) is a Langevin-like equation for a reaction-diffusion 
process where n is a coarse-grained particle density. Indeed, it has been ar- 
gued in [Sni that (fTH|) is the equation governing the pair contact process with 
diffusion, 2A —>■ 3A, 2A 0, the critical behavior of which is currently under 

^The difference between tlie Ito or tlie Stratonovich results is a trivial shift in /i. 
^Strictly speaking, at unbinding the system is in a state that evolves continuously to 
the absorbing state. 
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debate (see [S3] and references therein). A different realization of (|18|) comes 
from the field of synchronization in spatially extended systems. Pikovsky et 
al. inZl established that the difference field, ?7,(x, t), between two coupled- 
map lattices follows (fTH|) with c = 0. Two replicas of a coup led- map lattice 
are synchronized when n = 0, n 7^ corresponding to the asynchronous 
regime, and the active- absorbing, or the pinning-depinning, phase transition 
is analogous to a synchronization transition. As a last example, (|18|) has 
been studied in j2H] in the context of spatio-temporal intermittency. The 
unsuspected connections between these problems are illustrated in what fol- 
lows. 

The mean-field approximation to (|18|) consists of discretizing the Lapla- 
cian as l/2dJ2ji'^j where = n(xj, t) and the sum is over the nearest- 

neighbors of i. Substituting the values of the nearest-neighbors by the av- 
erage field, (n), a closed Fokker-Planck equation is obtained for P{n,t, (n)). 
This approach takes into account the effect of the noise and of a spatially 
varying order parameter. The steady-state solution is then found from the 
self-consistency requirement 

_ I^dn nPjn, (n)) 
^""^ S^dnP{n,{n))' ^''^ 

resulting in a phase diagram equivalent to that of figure (4) with a tricritical 
point at (/i = Q,hyj = 0) jHElEZI- Thus, for repulsive walls a complete wetting 
transition characterized by 2; = 2, i/^ = 1/2, and On = 1, is found at /Zc = 
for any 6 > 0. Different mean- field approaches yield different results for the 
exponent Pn, namely 1/p and D/2 ^35j, but this has been clarified in 
where a crossover from 1/p to a nonuniversal, continuous exponent D/2u 
was identified. For attractive walls, the three regimes previously described 
are recovered, namely, an active phase for /x < /i^ = 0, an absorbing one for 
/i > fi*{b), and phase coexistence of the active and absorbing phases in the 
region < /i < A qualitatively equivalent picture is obtained using 

a different space of parameters, for instance the strength of the noise vs. /x 

Beyond mean-field theory, (|18p is known to be superrenormalizable above 
6 = 0, i.e. the Feynman graphs may be computed to all orders and resummed 
IHl] ■ This does not imply that all the critical exponents are given, as 
the renormalization group-flow equation has runaway trajectories that are 
supposed to converge to the strong coupling fixed point. Representative 
results for d = 1 are summarized in table (2). 



18 



The rich phenomenology obtained for attractive walls is also observed 
within the MN framework. Active and absorbing phases coexist over a finite 
area /ic < A* < f^*{b), and lose their stability at, and fic, respectively. 

The phase-coexistence regime may be characterized by analyzing the single- 
site stationary density function, defined as the average of n{t) over pinned 
states rather than over all runs. This is depicted in figure (6) and shows how 
the histogram develops a maximum at n = as approaches the stability 
edge beyond which it changes abruptly into a delta function and the 

pinned phase becomes unstable. The simultaneous presence of two-peaks 
indicates that a fraction of the interface depins. This is clearly seen in the 
space-time snapshot of a numerical solution of (fTH|). on the left of figure 

(6) . The main features of such a pattern, spontaneous formation of domains 
with a wide range of sizes and lifetimes, has been identified as distinctive of 
spatio-temporal intermittency |58j . 

Numerical simulations provide evidence of a first-order phase transition at 
yUc- As was pointed out previously, both first-order and continuous transitions 
are observed at f^*{b) depending on the model parameters. For instance, a 
continuous DP-like transition was recently reported in [32] as the numerical 
solution of ()18p with b = — 9,c = 8, z/ = 0.1, and D = 1. By contrast, for 
b = —3, c = 2,1/ = D = 1, a first-order transition was found. The minimal 
Langevin equation that describes the DP universality class is 

^^1^' = z/V^ra - /in - bn"^ + ^Jnrj, (20) 

which differs from p8|) by the n^/^ noise amplitude. Thus, it appears that 
DP-like noise takes over the MN-noise, proportional to the field n, as the wall 
changes from repulsive to attractive in a given range of parameters. Deriving 
(j20|) from (jl8|) is an interesting theoretical problem that remains open. 

The similarity of behavior of KPZ interfaces bound by short- and long- 
ranged substrates also shows up in the framework of the order parameter n. 
Recall that an interface bound to the substrate is described by a particle- 
density field with low (high) density segments corresponding to detached 
(pinned) interface domains. Interface fluctuations can then be associated 
with the creation, annihilation, or merging of clusters of particles (see figure 

(7) ). In particular, the disappearance of triangles corresponds to a cluster- 
coarsening whose dynamics, being ultimately controlled by DP exponents, 
can be modeled by a contact-process. Recently, a lattice model of a gen- 
eralized contact process with long-ranged interactions between the edges of 
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low-density segments has been investigated jnij. A transition in the DP uni- 
versahty class is found for forces that decay sufficiently slow, and a first-order 
transition otherwise. Clearly, in terms of h this translates into a long-ranged 
interaction between the vertices forming the triangle bases, and it is reason- 
able to assume that, in turn, an effective long-ranged attraction between 
the substrate and the interface must obtain. Consequently, both short- and 
long-ranged interactions yield the same behavior below |51j . 

A direct application of the Cole-Hopf transformation n = exp{—h) to ()17p 
with A > in the presence of short-ranged interactions yields the Langevin 
equation 
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uV n - 2- an - hn - 2c7r + nri, (21) 

ot n 

which is (II 8|) plus the extra term — 2(Vn)^/n. Apart from the factor -2, 
equation (j2H) is the Cole-Hopf transform of 

dh 

— = V^/i + fi + be-'' + ce-^'* + r] (22) 

C/6 

which is the equilibrium EW model in the presence of a binding wall. Note 
that the factor of 2 in ()21|) cannot be absorbed by reparametrization. Indeed, 
equation PT|) does not admit proper mean-field solutions, it is not amenable 
to standard perturbative field theoretical methods, and cannot be integrated 
numerically [SB]. Nevertheless, (|2ip may be studied within the framework of 
active/absorbing phase transitions by simulating discrete growth models, and 
monitoring the order parameter n = exp{—h). This was done in |38j for two 
different models argued to correspond to a bound KPZ with A > and short- 
ranged interactions. Results suggest the existence of a new universality class 
characterized by /3„ = 0.33, ~ 0.22, z = 1.5,1^3, ~ 1. Ginelli et al. carried 
out a mean-field analysis based on a lattice model for KPZ-like growth, the 
single-step model, that does not completely disregard fiuctuations [HI]- They 
found a surface exponent 9n = 1/3 and succeeded in reproducing qualitatively 
the main features of the system for both signs of A. Finally, for a sufficiently 
attractive wall, the transition becomes first-order. 

h or nl The Cole-Hopf transform puts the problem of bound KPZ in- 
terfaces in a rather different context. For A < 0, a more compact equation, 
instability-free and amenable to a mean-field analysis, is obtained, whereas 
for A > the presence of (VnY/n suggests that the interface language is 
most natural. We stress, however, that the form of the continuum equation 
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does not determine the "correct" oder-parameter: In fact, proper scaling 
behavior is found for both h and n in discrete models. DP is a counterex- 
ample since the interface representation of the contact process shows clear 
signs of anomalous scaling Another interesting point concerns the n- 
representation of the EW problem, that is, the equilibrium limit (A = 0) of 
(jl7|) . The Cole-Hopf transform fixes the ratio \/v and thus the A = limit 
cannot be taken when the transformation is used. A comparison of published 
results on equilibrium and nonequilibrium wetting, reveals that the EW plus 
a wetting potential is not equivalent to the KPZ in the weak-coupling regime 
plus the same wetting potential. 

5 EXPERIMENTAL REALIZATIONS 

Since the idea of a wetting transition was introduced in 1977 by Cahn jM] , 
much experimental effort has gone into its realization, ultimately awarded 
by the confirmation of a large number of theoretical predictions. The ob- 
vious question then arises of what are the experimental realizations of the 
phenomenology described for the bound KPZ interface. As was discussed 
in section 2.2, the KPZ non-linearity models lateral growth. Although this 
mechanism is unlikely to be relevant in simple fiuids, it may be important in 
describing systems with anisotropic interactions for which the growth of tilted 
interfaces depends on their orientation. For instance, it has been shown that 
crystal growth from atomic beams when desorption is allowed is described 
by the KPZ equation [HH]. Furthermore, a renormalization group analysis re- 
veals that the KPZ term is always generated, except when excluded by sym- 
metry, whenever elastic objects depin in the presence of anisotropy [66j. This 
suggests that an answer may come from the field of crystal growth. Crystals 
grown by deposition of material onto a crystal substrate from a vapor phase 
can display wetting behavior. At sufficiently low temperatures, the adsorbate 
that intervenes between the substrate and the vapor may form small crystal- 
lites, which melt upon increasing the temperature, leading to the formation of 
liquid droplets. On approaching the wetting temperature, the liquid spreads 
over the substrate and coats it uniformly. Growth of solid phases is also 
possible, but in this case the wetting layer is likely to be under stress due to 
the misfit between the lattice constants of the substrate and the adsorbate. 
Two well-known relaxing mechanisms are the formation of misfit disloca- 
tions and the deformation of the surface (Asaro-Tiller-Grinfeld instability 
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[nZj). Faceting and, in particular, the appearance of pyramidal structures 
similar to those of figure (5) have been reported in simulations of wetting- 
layers and island formation in heteroepitaxial growth [HH]- These structures 
can either form directly on the substrate (Vollmer- Weber growth) or on top of 
a wetting-layer of finite thickness (Stranski-Krastanov-like growth). Clearly, 
further work is needed before a connection between crystal growth and a 
bound KPZ can be substantiated. 

A second direction of research is the experimental verification of syn- 
chronization transitions in extended, one-dimensional systems. The order 
parameter for such transitions is the difference between two dynamical sys- 
tems coupled to each other, n{x,t) = \ui{x,t) — U2{x,t)\. Starting from dif- 
ferent initial conditions, the systems become synchronized in the stationary 
limit, (n) = 0, above a certain critical coupling strength, and asynchronized, 
(n) > 0, otherwise. The character of the transition depends on whether the 
largest Lyapunov exponent of the system. A, is zero or negative. In the for- 
mer case, the transition is in the MN universality class, while in the latter it 
is DP-like or discontinuous depending on the system details. Equation (jl8|) 
is a stochastic model for the difference field n, with A = /ic at the transition. 
Examples of such chaotic coupled-systems include semiconductor laser arrays 
and liquid crystals describable by the anisotropic Ginzburg-Landau equation, 
for which the aforementioned behavior should be experimentally detectable 
[UnilZni- This goes hand in hand with finding a system which reproduces the 
critical exponents of directed percolation, a problem that so far has defied 
all attempts of solution despite extensive empirical efforts [7T] . 

The other large class of bound KPZ — short-ranged forces, positive non- 
linearity, and a lower wall — should also describe actual nonequilibrium wet- 
ting phenomena. There is nonetheless one other possible experimental real- 
ization. Sequence alignment is a powerful tool for determining relatedness 
between sequences of proteins or DNA segments. Algorithms have been 
devised to detect such correlations, the Smith- Waterman being among the 
most sensitive for finding related sequences in a database [72, . Scores are as- 
signed to each character-to-character comparison: positive for exact matches 
and negative if the two are different, or if pairing of an element with a gap 
occurs. An optimal alignment is one with the maximum possible score. Re- 
cently, Hwa et al. have mapped this algorithm to a KPZ with a lower wall, 
where the location of the critical point corresponds to the choice of the op- 
timal scoring parameters [75] . 
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TABLE CAPTIONS 



Table 1: Upper critical dimensions and critical exponents in the mean-field 
(MF) and the fluctuation regimes for complete wetting transitions with short 
and long-ranged forces. At the mean-field level, the exponents for short- 
ranged interactions are recovered in the m — > oo limit of the long-ranged 
ones. 

Table 2: Summary of the universality classes discussed in the text for short- 
ranged forces. For transitions other than DP, the theoretical predictions are 
z = 3/2, = 1, and /?„ > 1 if A < 0. 
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TABLE 1 



Exponent 


Short-range forces 


Long-ranged forces 




MF 


d<dc^2 


MF 


d < d^{m) = ^ 







1/4 


l/(2+m) 


(2-d)/4 




1/2 


2/3 


(2+m)/(2+2m) 


2/(d+2) 







1/3 


l/(l+m) 


(2-d)/(d+2) 




2 


2 


2 


2 
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TABLE 2 



Exponent 


A < 


A > 


DP ino] 


h 40^ 


n [55^ 


h 38 


n 38 


Oh,n 




1.1(1) 


0.355(15) 


0.215(15) 


0.1595 






1.0(1) 




0.97(5) 


1.7338 


Ph,n 


0.41 


1.5(1) 


0.52(2) 


0.32(2) 


0.2765 


Z 




1.52(3) 




1.55(5) 


1.5807 
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FIGURE CAPTIONS 



Fig. 1: Separation /i(x, t) of the AB interface from a rigid wall. 

Fig. 2: Typical effective interfacial potentials for negative and positive val- 
ues of b. Below the transition temperature, the potential well localizes the 
interface (dotted lines). 

Fig. 3: Left, pressure vs. temperature phase diagram of a pure substance. 
TP, CP and Ty^ stand for triple point, critical point, and wetting temper- 
ature; center, thickness of the wetting layer as a function of the relevant 
parameters for the paths indicated on the left and on the right; right, phase 
diagram in the b — fi plane. Solid and dotted lines correspond to continuous 
and first-order transitions, respectively (see text). 

Fig. 4: Phase diagrams for A < (left) from numerical solutions of ()17|1 in 
short-ranged potentials, and A > (right) from discrete models. The arrows 
denote different types of transitions explained in the text. 

Fig. 5: Snapshots of = 2, 1 interfacial configurations within the coexistence 
region, from numerical solutions of ()17|) in a short-ranged potential. 
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Fig. 6: Left, spatio-temporal evolution from a numerical solution of (|18|) 
within the coexistence region. Depinned regions {n < 1) are colored in dark 
gray and pinned ones [n > 1) in light gray. The instantaneous interfacial 
configuration for the time marked by the line is shown in figure (0). Right, 
single-site probability density function for different values of /i, also within 
the coexistence region. The stability edge is fi*{b = —4) ^1.3. 

Fig. 7: Typical interface profile within the coexistence region and its asso- 
ciated particle representation. A long-ranged attraction between the edges 
of low-density segments tantamounts to an a long-ranged attraction between 
the vertices forming the base of the triangle, which arguably induces an ef- 
fective long-ranged attraction between the substrate and the interface. 
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FIGURE 5 
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FIGURE 6 
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FIGURE 7 
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